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I. INTRODUCTION 

There is a rich literature on 6 He based on an effective a—n—n 3-body problem Besides 
pair interactions also ad hoc 3-body forces are used. The Pauli principle is approximately 
incorporated by projecting out "Pauli forbidden states" for the neutrons inside the a-particle 
wave function. While that approach catches presumably the halo structure of the two loosely 
bound neutrons clearly because of its strongly restricted ansatz it is not suited to probe 



ike the ones derived recently through 
21. 



modern nucleon-nucleon and three-nucleon forces, 
effective field theory and based on chiral symmetry 

Nevertheless some approaches already exist which directly attack the 6-nucleon prob- 
lem and beyond with realistic nuclear forces, namely in the method of no-core shell model 
(NCSM) 3 1 and the Greens function [4] Monte Carlo treatment. In [3j chiral two-nucleon 
and three-nucleon forces were used and applied to 7 Li with some under binding. In [4] the 
AVI 8 nucleon-nucleon interaction and a Urbana three-nucleon force was used again leading 
to some under binding, now for 6 He. There is also the stochastic variational Monte Carlo 
method jfl which, however, still applied simplified forces. 

The exploration of chiral forces goes on, also including explicitely the A-degree of freedom 
6J, which calls for an increased effort to establish rigorous approaches beyond A = 4. The 
achievments in [3] and [4] demonstrate that a direct treatment of 6 nucleons is feasible on 
present day computers and therefore we felt that another approach, the exact formulation 
within the Yakubovsky equations, is timely. About 20 years ago an analogous step turned 
out to be very fruitful, namely the exact formulation of the a-particle within the Yakubovsky 
scheme [7, 8]. This pioneering study opened the way to a nowadays standard treatment j^- 
121 ] and allows the inclusion of the most modern two- and three-nucleon forces and even first 
estimates of four-nucleon forces [3]. 

In section^ apply the Yakubovsky equations Q to the six-body problem using 
the basic notation for sub clusters [15| . In section HTT1 we add the identity of the nucleons 
which leads to a set of 5 coupled Yakubovsky equations related to 5 different sequential sub 
clusterings of 6 particles. In view of the expectation for the dominant structure of 6 He, 
namely an a-core and two loosely bound neutrons, we stop the sequential sub clustering 
with 3 fragments, though the additional step with two fragments could be easily performed. 

In section IIVI and the Appendices we provide technicalities which we consider useful for 
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a numerical performance. Finally we summarize in section |V] 



II. THE YAKUBOVSKY APPROACH TO 6 PARTICLES 

We use the standard notation a n to denote the various members of n fragments for a total 
of N particles. Here N=6 and a 2 , a 3 , a 4 , a 5 denote two-, three-, up to 5-body fragmentations. 
5-body fragmentations necessarily have one pair left and thus 05 can be used to point to a 
specific pair. 

03 C ct2 means that the three-body fragments 03 consist of sub clusters out of the two 
fragments in a 2 or a 3 D a 4 means that the sub clusters 03 when broken up lead to the sub 
clusters a 4 . 

The bound state \l/ obeys the homogeneous equation 

V = G J2Va 5 * (2-1) 

0.5 

where Go is the 6-particle free Greens operator and V a5 is a pair force. The first step is the 
summation of each pair force to infinite order. Defining 

^ a5 = G Q V as V (2.2) 

one obtains like for 3 particles 

i>a 5 = G t a5 ^ $a 5 b B 1pb 5 (2.3) 
bs 

where t a5 is a two-body t-operator obeying the Lippmann Schwinger equation 

ta 5 = V as + V a5 G t a5 (2.4) 



and S a5b5 = 1 - 5 a5b5 . 

Next one defines new components 

*Pa 5 a 4 = G t a5 Y ^aefcV'fcB ( 2 - 5 ) 
6 5 Ca 4 

where all pairs 05, b 5 are sub clusters of the fragments in a 4 . Clearly 

VWi = ^a 5 a 4 (2.6) 

a4D<25 



That relation (I2.6p is used to obtain a closed set of equations for t/j, 



0.50.4 ■ 



baC04 b4Dbs 

One separates now the components if) aBa4 for a given 04 from the rest : 

■0a 5 a 4 - G t a3 2j <5a 5 b 5 1pb 5 04 = G t as <5a 5 ft 5 E ^464^6564 ( 2 -8) 

&5Ca4 &5CQ4 64365 

Let us define for a fixed 04 the column vectors V^ 4 and -0( a4 ) with the components 

OTk = ^ 5 a 4 (2.9) 

and 

(^)) 6s = £ Wfcfc (2-10) 

b 4 Db 5 



Then introducing the matrix C ai with the elements C^ b5 = t a5 S a5 b 5 Eq. (12.81) reads 

(1 - G C a4 )^ a4 = G C a4 ifj ia4) (2.11) 

or 

^« = (1 - GoC a *)~ 1 G C f04 V' (a4) = G T a4 i/; M (2.12) 

Apparently T° 4 obeys 

T a 4 = C a 4 + C a4 G(jT a4 (2.13) 

In explicite notation (I2.12p and ( 12 . 1 3[) read 

^4 = Go £ T^f = G J2 T ak E (2-14) 

f>5Ca 4 6sCa4 &4D&5 

^565 = tarjasbs + E t ^5^o s c 5 G T^ b5 (2.15) 
csCa4 

Note, there are two types of T-matrices. For 0,4 of the type 123,4,5,6 T^ 4 b5 is a 3 x 3 
matrix and for of the type 12, 34; 5, 6 T® 4 b5 is a 2 x 2 matrix. 

Next we further decompose the right hand side of (I2.14p according to 3-body fragments 

a 3 : 

c = E E (2.i6) 

b 5 Ca4 b4Db 5 ,b4(Za 3 



and again 

0a 5 a 4 = Yl C (2- 17 ) 
0321014 

is an obvious consequence. 

Using again (I2.17P Eq. (12.161) can be rewritten as 

c = E G ° T «k E ^ E <Ih ( 2 - 18 ) 

Analogous to (I2.8P one separates b 3 = a 3 from b 3 ^ a 3 and gets 

c - G o E T Z, E w& = G ° E T ^ 5 E ( 2 - 19 ) 

&5Ca4 64Dfc5,fr4Ca3 &5CC14 043&5, &4Ca3 

where we defined 

vfe! = E wfe. ( 2 - 2 °) 

03=) 04 

Using a matrix notation it is easily seen that (12. 19[) leads to 

C = G o E E ^^Xt ( 2 - 21 ) 

b 4 Ca 3 6sC&4 

where D a a l a4 . hh obeys 

^os,a4;6 6 64 = ^565 ^4&4 + E E ^scJa^GoD^^.^ (2.22) 

C4Cd3 C5C04 

Note for a 3 = 1234, 5, 6 there are 18 pairs of a^, a±, for a 3 = 123, 45, 6 there are 9 pairs of 
05, (Z4 and for a 3 = 12, 34, 56 there are 6 pairs of a 5 , a^. This defines the dimensions of the 
different D-matrices. In the following the single nucleons in a n will no longer be displayed. 

III. IMPLEMENTATION OF THE IDENTITY OF THE NUCLEONS 

We start from (12.141) choosing the case as = 12 and 04 = 12, 34 and obtain 

/ _ n rnl2,3A 1 (12,34) n ^12,34 , (12,34) /Q x 

^12;12,34 — <J"0-t 12,12 ^12 + (jr O i 12,34 ^34 K 6 - 1 ) 

which according to (I2.10p is 

012;12,34 = GqT±2'i2 (012,123 + '012,124 + 012,125 + "012,126 
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+ V , 12;12,35 + "012;12,36 + ^12; 12,45 + 012;12,46 + "012;12,56) 
+ ^0^12,34 (034,134 + 034,234 + 034,345 + "034,346 

+ 034;34,15 + "034;34,16 + "034;34,25 + 034;34,26 + "034;34,56) (3-2) 

It is easily seen, going back to the definitions (12. 2p and (12. 5 p together with the antisym- 
metry requirement for the total state \I/ that 

034,134 + "^34,234 + "034,345 + "034,346 + "034;34,15 + "034;34,16 + 034;34,25 + "034;34,26 + "034;34,56 
= ^13-^24 (012,123 + 012,124 + 012,125 + 012,126 

+ 012;12,35 + 012;12,36 + 012;12,45 + 0>12;12,46 + 012;12,5e) (3.3) 



Therefore (I3.2p turns into 



"012;12,34 — ^o(^12,'l2 + ^12,34-^13-^24) (012,123 + 012,124 + 012,125 + "012,126 
+ 012;12,35 + 012;12,36 + 012;12,45 + 012;12,46 + 012;12,56) (3.4) 

The coupled equations (12.151) written out for a 4 = 12, 34 and a 5 or b 5 equal to 12 or 34 
yield when acting with P = P13P24 from both sides 

0^12,34 p _ ^12,34 
ri 12,12^ — 1 34,34 

p T 12,34 p _ T 12,34 (o r\ 

rl 12,34^ — 1 34,12 \ -°) 

Then defining 

^12,34 _ ^12,34 ^12,34 5 f0 R x 

T = r i2,i2 + ^12,34^ (3-6) 

it follows that T 12, obeys 

T 12 ' 34 = t 12 P + t 12 G PT 12 ' u (3.7) 
Therefore (13 .4p simplifies to 

012;12,34 = G T 12 ' 34 (^12,123 + "012,124 + 012,125 + 012,126 
+ "012;12,35 + 012;12,36 + 012;12,45 + 012;12,46 + 012;12,56) (3.8) 

Starting again from (12.141) but now for = 12 and = 123 one obtains 

./, _ r t 123 w/ 123) 4- r t 123 w, (123) 4- r r 123 70 123) 

^12,123 — 1*0-1 12,12^12 + ^O 1 12,23^23 + ^0-* 12,31 ^31 



= GqT 1212 ("012,124 + 012,125 + '012,126 

+ "012;12,34 + "012;12,35 + "012;12,36 + "012;12,45 + "012;12,46 + "012;12,56j 

+ ^0^12,23 (023,234 + 023,235 + 023,236 

+ "023;23,14 + "023;23,15 + 023;23,16 + 023;23,45 + "023;23,46 + 023;23,56) 

+ GoTi™ 31 (031,314 + 031,315 + 031,316 

+ 031;31,24 + "031;31,25 + 031;31,26 + 031;31,45 + "031;31,46 + 031;31,56) 

= Go(T^2 + ^12,23 P 12 -?23 + 7l2 2 ,31-Pl3-P23) ("012,124 + "012,125 + "012,126 

+ 012;12,34 + ^12:12,35 + 012;12,36 + 012;12,45 + 012;12,46 + 012;12,5e) (3.9) 

where we used permutation properties similar as in (13. 3p . 

The corresponding coupled sets (12.151) for a\ = 123 and using relations like 
-Pi3-P23^23-P23^i3 = *i2 reveals that 

t-i123 — 7^123 _i_ 7-1I23 r> JD 1 7^123 TD TD (o i n\ 

1 = 1 12,12 + 1 12, 23 ^12-^23 + 1 12,31 [6.W) 

obeys the equation 

T 123 = t l2 P + t l2 PG T 123 (3.11) 

where P = P 12 P 23 + Pi 3 P 23 - 
Then (13.91) simplifies to 

^12,123 = G T 123 (^12,124 + 012,125 + 012,126 
+ ^12;12,34 + ^12;12,35 + 012;12,36 + ^12;12,45 + ^12;12,46 + 012;12,5e) (3.12) 

The next step is to decompose 0i2,i23 according to ( 12.161) . For 05 = 12, 04 = 123 the 
possible a 3 's are: 1234 - 1235 - 1236 - 123, 45 - 123, 46 - 123, 56. 
Lets begin with 

"012^123 = GqTi 2 3 12 (012,124 + "012;12,34) 
+ ^0^12,23 ("023,234 + "023;23,14) + Co^l^ll (031,134 + 031;31,24) (3.13) 

Since 

"023,234 + "023;23,14 = -Pl2-P23("012,124 + 012;12,34) 

^31,134 + 031;31,24 = ^13^23 (^12,124 + V>12;12,34) (3.14) 
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(I3.13P simplifies according to (13.1 ip to 



= GoT 123 (^12,124 + ^12;12,34) (3.15) 

Similarily 

= G T 123 (^ 1 2,125 + ^12;12,35) 
^3 = Gor 123 (^ 12 ,126 + ^12;12,36) (3.16) 

Again using symmetry properties one gets 

= G T 123 ij 12 , 12A5 (3.17) 

^Sm = G T 123 iJ W6 (3.18) 

^i?i2 5 3 = G T 123 ^i 2il2> 56 (3.19) 

All summed up 

012,123 = VSs + + ^S'S + Vffig + l^g + itfS (3-20) 

agrees with (I3.12p . 

Next we decompose (I2.14p according to (12.161) for a*, = 12, = 12, 34. The possible g^'s 
are: 1234 - 125, 34 - 126, 34 - 12, 345 - 12, 346 - 12, 34, 56, which are now regarded in turn. 

;12,34 — ^0^12,12 ("012,123 + '012,124) + G T U ' M (034,234 + 034,134) (3-21) 

Since 

^13^24(^12,123 + 012,124) = 034,234 + ^34,134 (3.22) 

one can use (13. 6p and gets 

■0lf 3 12,34 = ^o(^12,'l2 + r^ 2 '!^ P13P24) (012,123 + "012,124) 

= GoT 12 ' 34 (^12,123 + 012,124) (3.23) 

Next 

/ 125,34 n ^12,34 , . n rpl2,34/ / . / \ 

12,34 — ^O 1 12,12^12,125 + ^O 1 12,34 (^34; 15,34 + 'V ; 34;25,34j 



s 



^12 ;12,34 — ^0^12,12 ('012;12,35 + ^12;12,45) + ^0^12,34 ^34,345 (3.24) 

The two amplitudes ^12-12 34 an d ^iS^ 34 can n °t be related by permutations, but their 
sum can be used 

,125,34 . ,345,12 n ^-,12,34/ , . , . , \ 

^12; 12,34 + fe; 12,34 ~ ^O 1 12,12 1^12,125 + ^12;12,35 + ^12;12,45j 

+ ^0^12,34 C034; 15,34 + "034;25,34 + ^34,345) (3.25) 

in the sense 

^13^24(^12,125 + ^12;12,35 + ^12512,45) = ^34;15,34 + ^34;25,34 + ^34,345 (3-26) 

This leads to 

^12; 5 12,34 + ^12; 5 12,34 = G T 12 ' U (^12,125 + ^>12;12,35 + V , 12;12,4s) (3.27) 

Similarily 

^12^4 + ^12,34 = C T 12 ' 34 (-012,126 + ^12;12,36 + V , 12;12,4e) (3.28) 

and finally 

, 12,34,56 _ n ^12,34 , . n ^12,34 , _ n ^12,34 , (n orA 

^12;12,34 _ 12,12 rl2;12,56 + ^O 1 12,34 V34;34,56 — ^O 1 ^12;12,56 \O.A)) 

Thus, Eqs.([3J23D, i^M) - <^M) , summarizes to 

, _ ,1234 , ,125,34 . ,345,12 . ,126,34 . ,346,12 . ,12,34,56 /, Qn \ 

^12;12,34 - ^12;12,34 + V ; 12;12,34 + ^12;12,34 + ^12;12,34 + V ; 12;12,34 + V ; 12;12,34 {0.6U) 

which when written out agrees with (13.81) . 

The two amplitudes V'lim an d ^12-12 34 expressed in ( 13. 15ft and ( 13.23[) are connected to 
each other as shown now. The expression (13.201) can easily be converted to ^12,124 and using 
in addition (13.301) one finds 



,,,1234 r rpVIZ, 1 1234 , / 1234 \ 
^12,123 _ ^O-f 1^12,124 + ^12;12,34J 

("i rpl2Z( 1 1245 , ,1246 , , 124,35 , , 124,36 . , 124,56 . ,125,34 
— trO 1 1^12,124 + rl2,124 + rl2,124 + ^12,124 + rl2,124 + ™12;12,34 

,345,12 , ,126,34 . ,346,12 . , 12,34,56\ fn oi \ 

' , 12;12,34 t 12;12,34 ' t 12;12,34 ' , 12;12,34/ [Q.Q1) 



Correspondingly ( 13.23j) yields 



? /,1234 _ ^ ^12,34/ ,1234 , , 
( / ; 12;12,34 ^O-f 1^12,123 + f 



1234 \ 
12, 124 J 
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GnT 12 ' 34 ^ 1235 + ^ 123 ' 45 + ^ 1236 + ^ 123 ' 46 + t/, 123 ' 56 + t/, 1245 
^(H I Y 12, 123 ^12,123 ' ^12,123 ' r 12,123 ' ^12,123 ' r 12,124 



/ 124,35 . , 
+ ^12,124 + ^ 



1246 , 1 .124,36 , 124,56 N 
12,124 ' ri2,124 "r ^12,124^ 



(3.32) 



with 



— P 34 , i/;J| 3 1 4 2 3 we can put (I3.3ip and f 1 3 . 3 2 j) into a matrix form: 



/ ?/; 1234 1 

^12,123 




|-G„( 


„/,1234 ■ 
\ r!2;12,34 / 





/ T 



123 



T 



12,34 



-T 123 P 34 T 123 



T 



12,34/ 



1 - P, 



,1245 , ,1246 
r!2,124 ' ri2,124 



31 J 



i / 124,35 , , 
+ ^12,124 + W 



,1234 
ri2,123 

,1234 
ri2;12,34 

124,36 , ,,124,56 , ,125,34 
12,124 ' ^12,124 ' ^12:12,34 



,345,12 ,126,34 ,346,12 ,12,34,56^ 

rl2;12,34 ' ^12:12,34 ' ^12:12,34 ' ri2;12,34 y 



/ 12;12,34 

, ,123,45 ,1236 , / .123,46 , ,123,56 , ,-, I" 
¥-12,123 ^ ri2,123 ri2,123 ^ ^12,123 ^ ^12,123 ^ ^12,124 

124,36 



,1235 
; 12,12: 

124,35 



123,46 
12,123 

124,56\ 



T ^12,124 T </"12,124 ^ ^12,124 ^ ^12,124 



(3.33) 



Since 



,,,1245 
ri2,124 

^124,35 
ri2,124 



■-^34^12,123 

TD I 123,45 
"-^34^12,123 



(3.34) 



the right hand side of (13. 33ft can be factored and f l3.33f) achieves the form 



/ ?/; 1234 1 

^12,123 






1 ? /,1234 ■ 
\ <Al2;12,34 / 





-T 123 P 34 T 



123 



1> 



1234 
12,123 



12,34/ 




,1234 
ri2;12,34 

123,45 



,1235 , ,12^,45 , ,1236 
^12,123 ' ^12,123 ' ^12,123 



i , 123,46 , , 123,56 
+ ri2,123 + ri2,123 

,125,34 . ,345,12 
rl2;12,34 i ri2;12,34 

i , 346,12 , 12,34,56 

\ + ri2;12,34 ~+~ ri2;12,34 



+ lb™' 3 * 
ri2;12,34 



(3.35) 



The right hand side can be reduced applying permutations and one obtains 

,1234 
¥1.2,123 

,1234 
y 12;12,34 



/ i/; 1234 1 

¥12,123 






i ?/; 1234 ; 

\ ¥12;12,34 / 





-T 123 P 34 T 123 





^34) 

-(P45 + iWg&3 + (1 " P 56 
(1 



D w , 125,34 
-H56 ) 1^12:12,34 



1 ,345,12 s 
+ ri2;12,34j 



123,45 
12,123 



1 , 12,34,56 
r "¥12:12,34 



(3.36) 



which can be put into the form 



,1234 
¥12,123 

,1234 
¥12:12,34 



with the result 



1234 
12,123 



•0 



1234 
12:12,34 




123,45 



v i"45 + ^46)^12123 + (1 — A6 — ^46)^12,12;! 



p w ,125,34 I , 345,12 
- r 56jl'¥l2;12,34 + ¥12;12,34 



)+1> 



12,34,56 
12:12,34 



(3.37) 
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where the matrix D obeys 




-T 123 A 4 T 123 



T 12,34 (l _ p \ 



34 
_ 34) 



-T lZ6 P u T 123 \ / D n D 12 , 

Go (3.38) 
T 12 * 34 (l - P34) J \D 21 D 22 

For the numerical treatment, however, we consider the structure (13.361) to be more ad- 
vantageous. 

The right hand side of (I3.36j) contains new amplitudes. After adequate permutation of 
(I3.30p one obtains from (I3.17P 

,123,45 _ s~i fT~?123 ( 1 1245 , ,123,45 . ,345,12 ,456,12 . ,126,45 . , 12,45,36\ /o oq\ 

ri2,123 ~~ (j i lri2;12,45 + ri2;12,45 ~+~ ri2;12,45 + ri2;12,45 + ri2;12,45 + ri2;12,45J (A^J 

or 

/ 123,45 rt rpYlZi /-i D \r> // 125,34 . ,345,12 \\ 
^12,123 - G T ("(I - P m) P Ml2-,i2,ZA + ^12;12,34)) 

= G T 123 (-P 3 5^2 2 ; 3 i 4 2,34 + tfggS) (3-40) 

Further (I3.27P yields inserting the decomposition of the right hand side related to (I3.20p 
and (E3nD yields 

,125,34 ^345,12 _ £^12,34/, 1235 , ^1245 , ^125,34 
rl2;12,34 ' ^12, 12,34 ~~ ^O- 1 1^12,125 ' rl2,125 ' ^12,125 

+ W.1256 ,125,36 ,125,46 ,1235 , ,124,35 , ,345,12 

^ ^12,125 ^ ^12,125 ^ r 12,125 ^ ( / ; 12;12,35 ^ ri2;12,35 ^ ( / ; 12;12,35 

+ W, 356 ' 12 + W, 126 ' 35 + W, 12 < 56 < 34 + ^ 1245 + t/, 123 ' 45 +7/; 345 ' 12 

( / ; 12;12,35 i ( / ; 12;12,35 ^ ri2;12,35 ' ri2;12,45 ' Vl2; 12,45 ^ ri2;12,45 

,456,12 ,126,45 . ,12,34,56 /q 

+ ri2;12,45 ' ri2;12,45 + ri2;12,45 l^- 4i J 

Here quite a few amplitudes can be related to previous ones by permutations leading to: 

,125,34 . ,345,12 ^ rr-,12,34 / r> D D D td \ I / 125,34 . ,345,12 \ 

fe; 12,34 + ^12,12,34 — ^O-* (-^35-^56 ~ "35 ~ -^46 ~ ^45 ) 1^12;12,34 + ^12; 12,34 J 

= GoT 12 ' 34 ((l — F 36 — P 34 — P46) ^34^35 ^12423 + (1 — -P34)-P46-P35' ? /'l2 3 123 

" (^45 + P35)^Tl 4 2,34 - PsstilflS) ( 3 - 42 ) 

Finally we regard V'12'12'34; which after (13.291) and suitable permutations of (I3.30p turns 
into 

,12,34,56 r T 12,34p p _/,12,34,56 
^12;12,34 ~ ^O-* i 35^46V , 12;12,34 
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= G T 12 ' 34 (P 3 5 J P46^ 1 1 g 4 2 ,34 + (1 - ^34)^35^46(^34 + < SI; ) ! (3-43) 

where we used 

^12;12,'56 = Co^ 12 ' 56 V ; 12;12,34 = PssPieGoP 12 ' 34 ^^;^^ = ^35^6^^2,34 (3.44) 

Thus we end up with 5 independent amplitudes: ^i 2 34 23 , ^i 2 34 2 ,3 4 5 ^i2, 3 23 5 t lhS,i2%i + 
^i2 5 i2,34> an d ^12; 12,34 > coupled in the equations 

-T 123 P 34 T 123 \ / i/j mii 



1 i/; 1234 1 

J ^12,123 






1 .1.1234. , 
\ ri2;12,34 / 





12,123 




p 34 ) / V ^ 



12:12,34 



c , ^ 23 (-P 34 ) T 123 \ / -(P 45 + P 46 )^ 234 3 + (1 " ^56 - . | . , 

T 42 > 34 (1 - P 34 ) J I (l-P^^ f> 



,123,45 ri rril23/ /-i d \D / , 125,34 . ,345,12 \\ 
^12,123 - (-(1 - ^36)^53(^1^,34 + ^12;12,34)) 

,1234 
; 12;12,34 



GoT 123 P 35 (^ 23 l 34 + 0ggg) (3.46) 



,125,34 .,345,12 rT 12,34/p p p p p V?/, 125 ' 34 J.,/, 345 ' 12 ^| 

^12;12,34 + ri2, 12,34 — ^O 1 {-^35^56 ~ -^35 — -<46 — -T45 ) IV ; 12;12,34 + ri2;12,34j 

G T 12 ' U ((1 — P 36 — P 34 — P 46 ) -P34-P35"0l 2 3 l 4 2 3 + (1 ~ -P34)-P46-f35' ! /'l 2 3 l 2 3 

(^45 + - P^i&SS) (3-47) 



,12,34,56 r T 12,34p p ,12,34,56 _ r ^12,34/ p p ,1234 
ri2;12,34 ~~ ^O-* - r 35- r 46^12;12,34 _ ^0^ \ r 35- r 46V ; 12;12,34 



+ (1 - ^34)^35^6(^2X34 + «2,3 4 )) (3-48) 

It remains to establish the expression for the total state We use ( 12. ip . ( 12. 5p . (I3.20p . 
(I3.30p together with permutations and obtain 



!;12,34j 



* = [1 - P23 - ^24 " Pl3 - Pu + ^13^24] [(1 " i^gua + ^gt 
— [1 — P 2 3 — P 24 — P13 — P14 + Pi3P 24 ] 

[((1 - P 34 )(P 45 + P 46 + P35 + Ae) - (^35^46 + PseP^M^s 

+ (-^45 + -P46 + -^35 + -P36 ~ -P35-P46)^l 2; 3 l 4 2 ,3 4 ] 

~~ [-P25 + -^26 + Pl5 + Pl6 ~ P13P25 ~ P13P26 — -^14-^25 — A4-P26 — -P15-P26] 
[(1-P3 4 )^S 4 23 + ^2 3 1 4 2,34 
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— ((1 - P34XA5 + Pas + P35 + Pm) - (PssPm + PsAM^s 

— (P45 + -^46 + -^35 + -^36 — -P35-Rl6)' ? /'l2; 3 l2,34] 

+ [1 - P 23 - P 24 " Pl3 - Pl4 + Pl3^24][(l - P34 - ^35 - A 6 )(l - Pm - Pl6)ti 2 2%l 

+ (1 — P45 — -^46 — -P35 — -P36 + -P35-P46)((l — -P56)(^12;12,34 + ' i /'l2;12,34) + ^12;12,34 )] 

— [P25 + P26 + Pl5 + Pl6 — P13P25 — P13P26 — -P14-P25 — -P14-P26 — P15P26] 

[(1 — -P34 — -P35 — -P36XI — P56 — P^^Pi^nl + (1 — -P45 — -P46 — -P35 — -P36 + -P35-R16) 

/ / -1 D w ,125,34 . ,345,12 \ . ,,12,34,56\i /o An\ 

((1 - F 5 6)(^12;12,34 + ^12;12,34) + ^12;12,34)] ( 3 - 49 ) 

The first piece 

#4 = [1 - P 23 " ^24 " Pl3 ~ PU + Pl3P24\ [(1 - J P34)^1 1 2 23 1 4 23 + VffiU] (3-50) 



Ala 



has exactly the form for a 4-nucleon bound state [7|, [9] based on the two Yakubovsky com- 
ponents %^i2^L23 an d '0i2 3 i2 34- Now for 6 nucleons these two amplitudes depend not only on 
the momenta of nucleons 1 to 4 but also on the momenta of nucleons 5 and 6. In the spirit 
of the effective 3-body model a — n — n ^4 factorizes into a product of the a-state and 
the wave function for the two neutrons. The next two pieces, still related to the same two 
Yakubovsky components, antisymmetrize the nucleons 5 and 6 in relation to the nucleons 
1 to 4. The remaining two pieces, going with the additional three Yakubovsky components 
^12,123; ^12-12 34 + ^12^2,34 an d '0i2;i2,34) §° clearly beyond the effective 3-body model and 
allow for additional sub clusterings of different types. In addition one has to keep in mind 
that the most general form (I3.49P allows many more distributions of neutrons and protons 
in the sub clusters than restricted in the effective a — n — n picture. 



IV. TECHNICALITIES FOR A NUMERICAL IMPLEMENTATION 



The 5 independent amplitudes require 5 different Jacobi momenta. For ^234 we choose 

at = ^(h- k 2 ) 

a 2 = ^(2A; 3 - h - k 2 ) 

a 3 = -fih -ki-%2- h) 

0*4 = -(£5 - h) 

a 5 = -(2(fc 5 + h) -h-k 2 -k 3 - fc 4 ) (4.1) 
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Then the individual momenta (under the condition Y2i k% = 0) i n terms of those Jacobi 
momenta are 



t - l - l - l - 
ki = at - -a 2 - -a 3 - -a 5 

k 2 = -ai - -a 2 - -a 3 - -a 5 

fa = a>2 ~ -a 3 - -a 5 

h = a 3 - ^a 5 

h = a 4 + ^a 5 (4.2) 



The kinetic energy is 

6 



EA_ = —(a 2 , + -al + -al + a\ + - ai) (4.3) 



2m m 4 3 

i=i 



The remaining choices of Jacobi momenta are given in Appendix |A] 
In a partial wave decomposition the basis states suitable for V'il 34 are 

\a1a2a3a4a5a1 > 

\(ti2-Mt 3 -)U(t 4 t 56 )TM T > (4.4) 

Here the orbital angular momenta k go with the two-body spins for particles ij, ji 

are total 1- and 2-body angular momenta coupled out of orbital and spin angular momenta, 
^3 and I4 are total 3- and 4-body angular momenta, I5 the total angular momentum of 
particles 5 and 6 relative to particles 1-4 and finally I4 and I5 are coupled to J, the conserved 
total 6-body angular momentum. The second state in f)4.4p refers to isospin in an obvious 
manner. 

Using (14. 4p the amplitude ip\ 23 4 has the representation 

l^iif 4 >= 5^ / n^ =1 daia*\aia 2 a 3 a4a5ai >< aia 2 a 3 a4a^a\\ip]^ > (4.5) 

where the discrete set of quantum numbers a\ runs over all values for a given J and T. 

Thereby antisymmetry requires that l\ + S12 + t\ 2 and U + S56 + ^56 have to be odd. 

Analogous basis states for the other 4 amplitudes and related Jacobi momenta should be 
obvious. One other example is given below. 
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If one projects the coupled equations (13.451) - f!3.48j) from the left onto the adequate basis 
states and expands the 5 amplitudes on the right hand side like in (14.51) one faces the task 
to evaluate the various kernels. As an example out of (13.451) we take 

< aia 2 a 3 a 4 a5Q!i|'0i2 3 i 4 23 > 

= -Go < a^a^^P^J^^M^^A > 

a l 

< a>24 a 4 a Xl^l2 3 123 > 

+ G < aia 2 a^a 4 a^ai\T 12Z \ ^ J ni =1 dbfy 2 \b[b 2 b 3 b'±b' 5 a' 2 > 



a 

,1234 



<b[b' 2 b^a' 2 ^ 1 \ M >+--- (4.6) 

where the • ■ • refer to the remaining operators and amplitudes and the b-states are related 
tn ?/; 1234 

lu Y 12-12 34 " 

Using techniques like the ones presented in [?], Q, l3, Q it is straightforward, 
though tedious, to generate the kernels like < aia^a^ailT^P^a^a^a^a^a^a^ > or 
< M 2 W 5 a 2 |T 12 ' 34 (l - P34XP45 + P46)|aia 2 a>4«5a / i >• 

In the course of the required recoupling among the different Jacobi momenta the variables 
for the 5 amplitudes on the right hand sides are in general linear combinations of intermediate 
integration variables which include angles besides momentum magnitudes. An example 
illustrates that situation: 

H = P 4 5^i 1 2 3 i23 (4-7) 
We project onto the basis states for Jacobi momenta of type b given in Appendix [A] 
\b> = \bib2b 3 b i b 5 ](l 1 s 1 2)ji{ks3A)j2(jij2)S(LS)I(ksse)j5(kj5)h(Ih)JM > 

|(tl2t34)tl-4(*l-4t56)TAf T > (4.8) 

where the orbital angular momenta k,i ^ 3 go with b{ and L goes with 63. The 2-body 
spins are Si 2 , S34 and s$q, I is the total angular momentum for particles 1-4 and I5 the total 
angular momentum of the pair 56 against the subsystem 1-4. They are coupled to the total 
angular momentum J, which is conserved. The isospin coupling should be obvious. 
For the sake of simplicity we choose s-waves which simplifies (14.81) to 

\b > = ^SI^s 12 ji^s S 4j 2 ^jsh^s 56 j 5 

\hhhhh > \(S 12 S 34 )S(SS 56 )JM > |(tl 2 t34)tl-4(tl-4*56)™ r > (4.9) 
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The basis states (14.41) related to the Jacobi momenta of type a and restricted to s-waves 

are 



\a> = 5 Sl2 j 1 5 j2 i8j A \8j hS s 6 8j h i 5 \a\a2aza±a 5 > 

\{si2-)h{h\)h{hsm)JM > \(t 12 -)t 3 (t 3 ^)U(t 4 t 56 )TM T > 



(4.10) 



Then ^i! 3 !^ has the representation 



| ^12,123 >= ' ' / nf =1 <iaja^|aia2fl304fl5 > \spin > a \isospin > a < al^ijjfm > (4-H) 

and H projected from the left is 



< b\H > 



< (tl2t34)*l-4(tl-4t56)rM r | < (s 12 S 3 4)5(5s 56 ) JM| < 6 1 6 2 b 3 6 4 fo 5 |P, 



E 



5 



f 11 
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I(^)^)*4(44)™t >< aWf 3 ^ > 
= 5 • ■ ■ S ■ ■ ■ I Ttf =1 daia%b < isospin\PH° spm \isospin > a 

b < spin\P^ m \spin > a < b^hb^P^^la^asa^ > 
< a|V>S*23 > 



(4.12) 



The permutation is separated into the 3 spaces: isospin, spin, and momentum, and the 
5's ■ ■ ■ are the strings of Kronecker symbols. 

The isospin- and spin-matrix elements can be calculated in a standard manner: 



b < isospin\P 45 p \isospin > a = 5 T T'hi T m' T K2t' 12 

^56 



\*ia+*i-4+l 



^3^34^1-4^56^4^56 



'12 



t' 



tl-4 ^34 



/ 1 

3 2 

1 1 

2 2 



t' t' 



T 



(4.13) 



6 < spin\P^ m \spin > a = S Sl2S ' i2 (-) Sl2+s+1 J P 3 s u Ss 5 6l'J 



I 5/ 
4 5 56 



8 1 V 

Sl2 2 J 3 
2 'S' 534 



J 3 2 



5 



2 2 S56 



j/ „/ 



4 S 56 



J 



(4.14) 
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For the momentum space part we insert a complete basis and obtain 

< hbzhhhlP™ ™ = J d% ■ ■ ■ d% < hhhhhlbfyb&b's >< b'^b'M'^ 

= (-L) 5 / db[ ■ ■ ■ d% < b{b\b^U^b\b^\P™ m (4.15) 
V4vr J 



The transposition P^ om acting to the left leads to a state with the same quantum numbers 
but different meaning. Particles 4 and 5 are interchanged. That state can be reexpressed 
again in terms of the old b-state using the relation between the Jacobi momenta &i, • • - 65 
and the Jacobi momenta with particles 4 and 5 interchanged. It results 

^^^^^ 1 1^ - 3^ 1 3^ - 3^ 

< &1M3M5I-P45 = < 61, -(62 - -63 - 64 - -65), + - 64 - -65), 

~ \h + 64 - Ik), -(b 2 + l -h + 64 - (4.16) 

Consequently 



< hb 2 b 3 b 4 b 5 \PS om = (^U 5 / A ■ ■ ■ db' 5 

V47T J 



a 1 1 - a 3-1 - 3- - 3 a 

< 616;, -(6 2 6 2 - 2 6 3 & 3 - b A - ^6565), -(-6 2 6 2 + -6363 - b 4 b' 4 - -hb' 5 ), 

l -{-b 2 v 2 - h 3 y 3 + - lhb' B ), -ib 2 v 2 + h 3 y 3 + h% - ±b&)\ (4.17) 

The state \b\ ■ ■ • b§ > can be reexpressed in terms of the state |ai ■ ■ • a 5 > as 

\hhhhh >= \h, ^(fe 2 - h), ~{2b 2 + b 3 ),h, b 5 > (4.18) 

Consequently using (14. 17ft and ( 14. 18ft one obtains 

< b\b 2 b 3 b^b 5 1 P 4 ™ om J TL^ =1 daiaf \ a\a 2 a 3 a^ >< aia 2 a 3 a^a$\^]^i 23 > 

= jdb[--- d% < b u ^\2b 2 h' 2 - 2b 3 %\, l -\ l -b 2 y 2 + - A b 3 y 3 - - h 5 k 
l\-\b 2 y 2 -h 3 y 3 + b 4 y 4 -h 5 y 5 )\, 

\b - 2b' 2 + h 3 b' 3 + hb' 4 - \hb'MSi23 > (4^) 5 (4-19) 



5J\i 



The angular integration over b[ yields directly An. Further one can put b 2 into the z- 
direction and b' 3 into the x-z plane. Therefore 

< bib 2 b 3 bib 5 \P™ m J U^ =1 daia^\aia 2 a 3 a4a 5 >< aia 2 a 3 a A a^\ilj^i 23 > 

= v^y dcoso' 3 db' 4 dy 5 < h, l -\2b 2 y 2 - 2b 3 y 3 \^\ l -b 2 y 2 + h> 3 y 3 - h A y 4 - h 5 y 5 )\, 
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\\ - \hV 2 - \b& + hb' A - l&A + i&A + - \hV,U%%* > (4.20) 

This is an example where a 4-dimensional interpolation appears necessary. Therefore if 
the momentum magnitudes are discretized choosing for each one a certain grid, interpolations 
are inevitable. 

In the 3- and 4-nucleon problems cubic Hermitean spline interpolation turned out to be 



very efficient 171 ] . Thus a 5-dimensional interpolation, for instance, has the form 



/r/fc2m(Ol<22&3 a 4«5) — 
3 3 3 3 3 

^ ^ 5 , r (ai)S , s (a 2 )S' t (a3)S' u (a4)S'„(a5)/(ai T .a2 S a3 t a4 ll a 5 „) (4.21) 

r=0 s=0 i=0 u=0 u=0 

where ijklm denotes the 5-dimensional cubus around the point G.1Q.2Q.3 040,5 and a ifc denotes 
the grid points for the variable 04. 

For ax, for instance, the 4 grid points related to are aio < (in < «a < 0-12 < 0-13 and 
similar for the other variables. Here for the sake of a simpler notation we renumbered the 
grid points in that context. 

Further / is the function to be interpolated and fijkim the interpolating one. The spline 
functions are given in Appendix [B] and the conditions underlying that form ( 14.211) can be 
found in 

The coupled set ( 13 .45H - ( T3T48T) in a matrix notation has the schematic structure 

r](E)^ = K(E)ip (4.22) 

where E is the searched for energy eigenvalue at which the auxiliary kernel eigenvalue rj(E) = 
or 3- and 4-nucleon bound states a Lanczos type algorithm turned out to be very efficient 



12 



181 ] . Starting from an arbitrary initial x/j — ifio one generates by consecutive applications 
of K a sequence of amplitudes ip n , which after orthogonalisation form a basis into which ijj 
is expanded. It turn out that a reasonably small number of K-applications (of the order of 
10-20) is sufficient, which leads to an algebraic eigenvalue problem of rather low dimension. 
Then the energy is varied such that one reaches t](E) = 1. 

If one regards the sub clustering underlying the 5 amplitudes only two of them, ^f 3 ^ 
and V'i2- 3 i2 34' are related to the very approximative effective 3-body model of an inert a-core 
and two neutrons. The total isospin quantum numbers of the 6 He ground state are T = 1 
and Mt = — 1 (if we define the magnetic isospin quantum number of the neutron as —\)- 
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Now even these two amplitudes depending on the Jacobi momenta of type a and b, (14. ip 
and flAip . respectively, are not restricted to £4 = and £56 = 1 in case of V'iim but also 
allow £4 = 1 and t 5S = 1, £ 4 = 2 and = 1 and £4 = 1 and = 0. Similarily the 
" deuteron-deuteron" like substructures of "0i|. 3 i2 34 are n °t restricted to £12 = £34 = and 
£56 = 1 but also other two-body isospins are allowed which do not built up a £ = en- 
core. The amplitudes ^12^23 g°i n g with Jacobi momenta of type c, (1A4[) . refer to a 3-body 
together with a two-body sub clustering, which is also not present in the effective 3-body 
model. The linear combinations tp^ 5 ^ 34 + ih^vz 34 refer again to a 3-body together with 
a 2-body sub clustering, where the underlying fragmentation related to 2-body fragments 
differs from V'lin.W And again this is beyond the effective 3-body model. Finally tpl^'n'sA 
allows for several additional 2-body sub clusters which are not contained in the effective 
3-body model, either. 

If one would add another step in the Yakubovsky scheme, namely to 2-body fragmenta- 
tions (Z2, the resulting amplitudes if)^ ai a3 would point to 5-body subclusters together with a 
single nucleon or to two 3-body subclusters (like 3 if — 3 H) in addition to 4-body and 2-body 
subclusters. Now all those additional structures are of course also generated by the coupled 
system (13 .451) - (13 .481) ending with 03 fragmentation, which we presented. In other words that 
Yakubovsky scheme is complete and delivers an exact description of the 6-nucleon problem. 



V. SUMMARY 



The Yakubovsky equations have been derived long time ago by O. A. Yakubovsky [14]. 
Therefore the application presented here could have been given also long time ago. It is, 
however, only now after the experiences with 3- and 4-nucleon problems in the Faddeev- 
Yakubovsky schemes that the technical expertise has been developed in the last decades 
and the very strong increase of computer power just recently achieved allows to attack 
the 6-body problem in that exact formulation. Therefore we felt it is timely to work out 
that scheme for that system. Another argument is the development of nuclear forces in 
a systematic manner in the realm of effective field theory and based on chiral symmetry. 
Two-, three-, and four-nucleon forces have been derived consistently to each other and they 
are waiting to be applied in light nuclear systems and checked against nature. Several tests 

nnnn 

in that spirit already appeared [3|, |4j, [12|, [19] but for the purpose of benchmarking the exact 
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approach in the Yakubovsky scheme is strongly recommended. 

Here we restricted the formulation to two-nucleon forces only but the inclusion of three- 

n 

nucleon forces can easily be done like pioneered in |8|. 



In section [TT] we used the general basic and standard formulation, where a n points to 
n-body fragmentations for a system of N > n particles. 

We worked out that scheme ending with in the spirit of the usually applied though 
approximate effective a — n — n 3-body model. 

Though the step to 02 could be easily done, ending with still includes exactly the 
whole dynamics. 

The Pauli principle is then exactly incorporated in section Uni leading to a set of 5 coupled 
equations for 5 independent Yakubovsky components, which built up the total state. 

The technical performance in a partial wave decomposition is only touched in section 
IIVI The 5 different Jacobi momenta as well as a necessary mult i- dimensional interpolation 
scheme, like modified cubic Hermitean splines, are given. For solving the high dimensional 
energy eigenvalue problem of the 5 coupled equations we point to the Lanczos type algorithm, 
which turned out to be very efficient in the 3- and 4-nucleon problem. It remains to work 
out the partial wave projected kernels, which is straightforward and can be carried through 
along the lines cited above. An example for that is presented. 

We expect that on the most modern supercomputers with parallel architecture this for- 
mulation can be numerically mastered. 
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Appendix A: Jacobi momenta related to the independent Yakubovsky components 

Here we display various Jacobi momenta related to the independent Yakubovsky compo- 
nents. 
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To ^12:12,34 belongs 

1 -> 

h = -(fa - fa) = oi 
1 - 

^2 = - fo) 

^3 = 7}(fa + fa~fa- fa) 

fa = ^(fa - fa) = a 4 

65 = -(2(fa + fc 6 ) - fa - fa - fa - fa) = a 5 
The individual momenta are expressed in terms of those Jacobi momenta: 





1-. 




fa 


= 6l + -&3 










— * 

fa 


1- 




= - 6 ' + 2' 
















fa 


= b 2 b-i 








4 5 








fa 


= - i« 










fa 


= &4 + 





The kinetic energy is 
To ^i2 3 i23 belongs 

Ci = ^(fci - fc 2 ) = «i 

c 2 = ^(2fc 3 - fa - fa) = a 2 

c*3 = fa) 
1 

c 4 = -(fc 4 + A; 5 - 2fc 6 ) 

C5 = \(fa + fa + fa - fa - fa - fa) 



fa 
fa 
fa 



Cl - -C2 - -C 5 
1 - 1 - 



-Cl 



2 C2 " 3 C5 



c 2 - 305 
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«4 = C 3 + -C 4 + -C 5 

r ^ U U 

«5 = -C3 + + g C 5 

A; 6 = -C4 + -c 5 

and the kinetic energy is 

V^. - — r 2 + — r 2 + — r 2 + — r ; 

' 2m m 4m m 4m 

i=i 

To V^fi^ belongs 

1 ,-> 

di = -(fci - Afc) = ai 
1 - 

d 2 = -(2fc 5 - fci - ^2) 

4 = ^(£3 - fc 4 ) = b 2 

d A = ^(k 3 + h- 2A; 6 ) 

db = ^(h + k 4 + k 6 -ki-k 2 - 



or 



1 - 1 - 

fci = di - -d 2 - -d 5 

. 1. 1. 
k 2 = -di - -d 2 - -d 5 

7 7 1 t 1 r 
h = d 3 + -rf 4 + -d 5 

1 1 - 

& 4 = -d 3 + -d 4 + -d 5 

h = d 2 - -d 5 

7 7 1 T 

k 6 = -d 4 + -d 5 



and the kinetic energy is 

y2^ = -di+^di+-di+^d 

' 2m m 4m m 4m 

To ^S'2 2 34 belongs 

ei = -(fc 3 - £4) = ^2 

e 2 = -(2fc 5 - fc 3 - ^4) 
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or 



e 3 = ^(fci - k 2 ) = Si 
e 4 = -(h + £ 2 - 2^) 

e 5 = ^(fci + k 2 + h-k 3 -k 4 - k 5 ) (AlO) 



«i = e 3 + -e 4 + -e 5 

^2 = -e 3 + -e 4 + -e 5 

n - i- i- 

^3 = ei - -e 2 - -e 5 

k 4 = -ei - -e 2 - -e 5 

«5 = e 2 - -e 5 

4 = -e 4 + ^e 5 (All) 



and the kinetic energy is 



E ^ = -c? + ±4 + lei + ^ + J-e§ (A12) 
2m m 4m m 4m 3m 

i=i 



rp , 12,34,56 i i 

To ^i2;i2,34 belongs 



k 2 ) 


= ai 


h) 


= 6 2 


fee) 


= a 4 



or 



h = \{h 

h = \{h 
1 

h = -{h + k 2 - k 3 - k A ) 

h = ^(2(4 + h) - & - fc 2 - fc 3 - fc 4 ) = & 5 (A13) 



fci = fi + \h~\h 
h = -fi + \h~ \h 

h = f 2 — ~U — -^h 

- 1 - 1 - 

k 4 = —fi - -J4 - -^h 

h = h + \h 
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(A14) 



and the kinetic energy is 



6 b 2 
2m 



i=i 



—fi + -ft + -fi + ^~fl + ^f! 

m m m 4m 8m 



(A15) 



Appendix B: Modified spline functions 



Choosing four grid points x , X\, x 2 , and x 3 such that x\ < x < x 2 the modified spline 



functions 



17] are 



with 



S (x) 
Si(x) 
S 2 (x) 
S 3 (x) 



x 2 — X\ 



-fa{x 

X\ — Xq X 2 — Xq 

, . i I \ I "^l %1 -^0 \ 



MX 



x 3 - x 2 



X\ — Xq X 2 — X\ X 2 — Xq X 2 — X\ X3 — X\ 

Xi — Xo 1 , , ,X 3 — X 2 X 2 — X\ 1 

<fo{x) + h{x)z —~ — + M x >{- 



Mx 



X 2 — X\ X 2 — Xq 
X 2 — X\ 1 



X 2 — X\ X3 — X 2 X3 — X\ 



x 3 - X 2 X 3 - X\ 



0i (ac) 



4> 2 (x) 



4> 3 (x) 

(j) 4 (x) 



(x 2 - x) 2 

x 2 - Xx) 3 
(xi — x) 2 



x 2 - Xi) 3 

(x — xi)(x — x 2 ) 2 

(x 2 - Xx) 2 
(x — Xi) 2 (x — x 2 ) 

(x 2 - Xx) 2 



((x 2 - x x ) + 2(x - xi)) 
((x 2 - xx) + 2(x 2 - x)) 



(Bl) 



(B2) 
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